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We investigate the causal temperature profiles in a recent model of a radiating star 
undergoing dissipative gravitational collapse without the formation of an horizon. It is 
shown that this simple exact model provides a physically reasonable behaviour for the 
temperature profile within the framework of extended irreversible thermodynamics. 
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1. Introduction 

The Cosmic Censorship Conjecture has continued to occupy center stage within the 
realms of relativistic astrophysics. The final outcome of the gravitational collapse of 
a star is still very much open to debate with the discovery of models admitting naked 
singularities!"^ Various scenarios of gravitational collapse have been considered in 
which the energy momentum tensor is taken to be e ither a perfect fluid or an 
imperfect fluid with heat flux and anisotropic pressureP^El jtj is well known that 
the collapse of reasonable matter distributions always lead to the formation of a 
black hole in the absence of shear or in the case of homogeneous densities. It has 
been shown that shearing effects delay the formation of the apparent horizon by 
making the final stages of collapse incoherent thus leading to the generation of 
naked singularities.® In this letter we revisit a radiating stellar model proposed by 
Banerjee et aP, (hereafter referred to as the BCD model) in which the horizon is 
never encountered. The interior matter distribution is that of an imperfect fluid 
with heat flux and the exterior spacetime is described by the radiating Vaidya 
metric.^ The junction conditions required for the smooth matching of the interior 
and exterior spacetimes across a four-dimensional time-like hypersurface are solved 
exactly. 
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In this letter we investigate the physical viability of the BCD model. In par- 
ticular, we analyse the relaxational effects on the temperature profiles within the 
framework of extended irreversible thermodynamics. We are in a position to obtain 
exact solutions to the causal heat transport equation for both the special case of 
constant collision time as well as variable collision time. Our results are in agree- 
ment with earlier thermodynamical investigations of radiating stellar models. We 
find that relaxational effects enhance the temperature at each interior point of the 
stellar configuration. Our investigations show that the BCD model displays physi- 
cally reasonable temperature profiles throughout the evolution of the star. 



2. The BCD radiating model revisited 



In the BCD model the following form of the metric for the interior spacetime is 
assumed 



ds 2 = -A 2 (r,t)dt 2 + B 2 (r,t)[dr 2 + r 2 d9 2 + r 2 sin 2 



(1) 



in which the metric functions A and B are yet to be determined. The energy mo- 
mentum tensor for the interior matter distribution is given by 



r"" = {p + P) vV + pg^ + gV + q v v>* 



(2) 



The heat flow vector g M is orthogonal to the velocity vector so that q^v^ = 0. In 
order to generate an exact model of radiative gravitational collapse, the following 
ansatz was adopted for the metric functions inQ 

A = o(r) (3) 
B = b(r)R(t) (4) 

which reduces the Einstein field equations for the interior matter distribution to 
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where 'dot' and 'dash' indicate derivatives with respect to time and the radial 
coordinate respectively. The condition of pressure isotropy yields 
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(8) 



Since the star is radiating energy the exterior spacetime is described by the Vaidya 
metric given explicitly in the form 

2M(v) 



ds' 



1 



dv 2 - 2dfdv + f A {dB 2 + sin 2 



(9) 
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where v is the retarded time and M(v) is the exterior Vaidya mass. The junction 
conditions required for the smooth matching of the interior metric Q and the 
exterior Vaidya metric © across a time-like hypersurface S are given by 



(?\B)s = r s 

r 3 BB 2 
2A 2 



r 2 5' 



25 



(10) 
(11) 

(12) 



where represents the total mass of the stellar configuration of radius r inside 
E. Utilising @ and 10 in the boundary condition 1111(1 yields 



2RR + R 2 + mR = n, 
where m and n are constants. A simple particular solution of (|13|l is 

R(t) = -Ct, 



(13) 



(14) 



where C > is a constant of integration. As pointed out in Ref0the mass-to-radius 
ratio, ms/fs, is independent of time. A simple calculation yields 
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(15) 



where b(ro) = 6o and ro dehnes the boundary of the stellar configuration. It is 
interesting to note that the parameters in (|15fl may be chosen so that 2ms /fs < 1 
in order to avoid the appearance of horizon at the boundary. 



3. Causal temperature profiles 

In this section we consider the physical viability of the BCD model. In order to 
satisfy the condition of pressure isotropy ©, the BCD model assumes b(r) = 1 and 

A = a(r) = (1 + ^or 2 ) (16) 
which in view of (17-20) leads to 

m = -4£ r , n = 4£o(l + £ r 2 ), C = ~ ( - \m\ + (m 2 + An) 1 ' 2 ^ . (17) 

The fluid volume collapse rate is 

e ~AB-(l+t a r 2 )t (18) 

which is the same in both the radial and tangential directions in the absence of 
shear. The proper stellar radius is given by 

r p (t) = / Bdr = -Ctb . (19) 
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Since the star is collapsing we require that C be positive which corresponds to 
— oo < t < 0. We further have 



C 2 < 4£ (1 +£or 2 ). 
The Einstein field equations ©-10 reduce to 

3 



t 2 (l+Zor 2 ) 2 
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t 2 (l+£ a r 2 ) 2 
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(22) 
(23) 



(1 + Cor 2 ) 2 CH 3 

We note that all the above thermodynamical quantities diverge as t — * 0. The 
regularity conditions p > 0,p > and p' < 0, and p' < together with the dominant 
energy condition, (p — p) > and the more stringent requirement (p + p) > 2\q\ are 
all satisfied when 

2^ 2 



1 - 
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We can now write 



We note that that when 
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C 2 r 2 
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C 2 <l/r 2 +£ 2 r 2 + 2£ 



(24) 



(25) 



(26) 



the boundary surface can never reach the horizon. El Furthermore, the surface red- 
shift is given by 

(27) 



1 + *jj= ( l + ro-^+ro&o 



which diverges for an observer at infinity at the time of the appearance of the 
horizon. For the BCD model Q27JI reduces to 



i + z s = (i - cvor 1 



(28) 



which diverges when C = l/rg. In order to avoid the divergence of the surface 
redshift we must have 

l/r 2 < C 2 < l/rl + e r 2 + 2£ , (29) 

where we have taken (126(1 into account. The luminosity of the star as perceived by 
an observer at infinity is given by 

dm 
dv 



L = 



c 3 r 3 (l + £ r 2 -rC) 
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(30) 
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which is independent of time. We now turn our attention to the evolution of the 
temperature profiles of the BCD model. To this end we employ the causal transport 
equation for the heat flux, which in the absence of rotation and viscous stress is 
given by 

Th a b q b + q a = - K (D a T + Tu a ), (31) 

where r is the relaxation time for the thermal signals. Setting r = in the above, we 
regain the so-called Eckart transport equations which predict infinite propagation 
velocities for the dissipative fluxes. For the line element the causal transport 
equation l|31|) reduces to 

T( q B), t + A(qB) = -K < ^^ (32) 

which governs the behaviour of the temperature. Setting r — in il-il'l we obtain 
the familiar Fourier heat transport equation 

A(qB) = - K ^|k (33 ) 

which predicts reasonable temperatures when the fluid is close to quasi-stationary 
equilibrium. In order to study the evolution of the temperature in the BCD model 
we employ the thermodynamic coefficients for radiative transfer as outlined in Ref 9 
The thermal conductivity takes the form 

k = 7 T 3 r c , (34) 

where 7 (> 0) is a constant and t c is the mean collision time between the massless 
and massive particles. We further adopt the generalised power-law behaviour for r c 



r c = T" w , (35) 

where a (> 0) and gj (> 0) are constants. The velocity of thermal dissipative signals 
is assumed to be comparable to the adiabatic sound speed, which is satisfied if the 
relaxation time is proportional to the collision time: 



(36) 



where r (> 0)is a constant. The constant (3 is a measure of the strength of relax- 
ational effects, with (3 — giving the noncausal case. Using the above definitions 
for r and k, (|32|l takes the form 

P{qB),tT~ u + A(qB) = _ a T3 -"(AT),r _ (37) 

The Eckart temperature is readily obtained by setting /3 = in (|37|) . We are in a 
position to integrate (|37|l for the special case u> — which corresponds to constant 
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collision time and more interestingly, the case u> = 4 which gives a variable collision 
time. For constant collision time (lu = 0), the causal temperature profile is given by 



T\r,t) 



8/3^0 [2(rg - r 2 ) + (o(4 r 4 )] 8£„ [3(r 2 - r 2 ) + 3£ (r 4 - r 4 ) + £ 2 (r 6 - r 6 )] 
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where L is given by H30|l and S is a constant. For u> = 4 , the causal temperature is 
given by 

T 4 (r, i) 
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where the exponential integral Ei(z) is defined as 

poo e -t 

Ei(z) = - / dt. 



(40) 



We note that the noncausal temperature (/3 = 0) and causal temperature are equal 
at the boundary (r = vq). Figure 1 shows that the relaxational effects are dominant 
when the stellar fluid is far from equilibrium (large values of (3). In the case of 
variable collision time, figure 2, we see that the causal temperature is everywhere 
greater than the corresponding noncausal temperature within the stellar interior. 
Furthermore, figures 1 and 2 indicate that the causal temperatures at late times 
(large values of j3) decrease more rapidly than the causal temperatures when the 
star is close to quasi-static equilibrium. This is in agreement with the perturbative 
results of Ref PH as well as the acceleration- free model studied in Ref. 11 



(38) 



. + at(l + £ r 2 )] 



(39) 



4. Concluding remarks 

We have investigated the physical viability of the BCD model within the framework 
of extended irreversible thermodynamics. We have shown that this simple model 
allows us greater insight into the evolution of the temperature for different collision 
times. More importantly, we were able to confirm earlier findings that the causal 
temperature dominates the Eckart temperature within the stellar core, even for 
variable collision time. As pointed out in earlier treatments, the constant collision 
time approximation is only valid for a limited period of the stellar evolution!^ One 
expects that the collision time between the particles making up the stellar fluid to 
change with temperature. Such effects on the evolution of the temperature profiles 
were clearly demonstrated with the variable collision time solution. It must be 
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Fig. 1. Temperature profiles for constant collision time, (close to equilibrium - solid line), (far 
from equilibrium - dashed line) versus r. 
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Fig. 2. Temperature profiles for variable collision time, (close to equilibrium - solid line), (far 
from equilibrium - dashed line) versus r. 



pointed out that the truncation of the transport equations leads naturally to an 
implicitly defined temperature law. Such a temperature law may only be valid for 
a limited period of collapse. What remains is to investigate the behaviour of the 
temperature by employing the full transport equation for the heat flux as well as 
to include the effects of shear. The general framework for such an investigation has 
recently been provided in Ref ll2l for a spherically symmetric radiating star. 
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